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We study properties of phase transitions of 2D superconductor liquid crystal phases, and analyze 
the competition between the recently proposed Pair Density Wave (PDW) and nematic 4e super- 
conductor (4eSC). Nematic fluctuations enhance the 4eSC and suppress the PDW phase. In the 
absence of lattice effects, the PDW state exists only at T = and the low temperature phase is a 
nematic 4eSC with short ranged PDW order. A geometric description of the 4e SC is presented. 
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A beautiful series of experiments [l|,yl have shovifn that 
the cuprate superconductor La2-xBaa;Cu04 exhibits a 
remarkable dynamical layer decoupling behavior near the 
X — 1/8 "anomaly". In this regime the onset of static 
charge and spin stripe order coincides with the develop- 
ment of an extreme transport anisotropy. Similar effects 
have been seen in stripe-ordered Lai.6-2:Ndo.4Sr^Cu04 
[3, |j| and in the magnetic-field induced stripe-ordered 
phase of La2-2:Sr2:Cu04 5]. 

The remarkable layer-decoupling effect suggests that 
the inter-layer Josephson coupling is somehow frustrated 
when SC order is forced to coexist with charge and/or 

tin stripe order. It was proposed by Berg et al [Q- 
, that this effect can be understood naturally by pos- 
tulating that the SC order also becomes "striped" and 
that all three orders rather than competing are inter- 
twined. In the resulting striped SC state, a pair-density- 
wave (PDW) , the (unidirectional) SC pair field oscillates 
in space with an ordering wave vector Q. The PDW 
has the same structure as an "Fulde, Ferrell, Larkin, 
Ovchinnikov" state, (FFLO)[^, [13) except that in this 
case the spatial modulation of the Cooper pairing is not 
due to Fermi surfaces mismatched by a (Zeeman) mag- 
netic field (as in the standard FFLO scenario). Instead, 
the PDW arises in a strong coupling regime where the 
BCS mechanism is not effective: the modulation of the 
SC state is due to the same physics behind the formation 
of inhomogeneous states in doped Mott insulators [ 111 [ia | . 
FFLO states have been proposed to occur in imbalanced 
cold atom Fermi systems with different species [13| and in 
heavy fermions systems when different orbitals hybridize 
under external pressure [lj|. In Ref. [l^ the phase di- 
agram of this problem was studied deep in the (Ising) 
nematic phase assuming that the coupling to the lattice 
is so strong that it completely suppresses the fiuctuations 
of nematic (orientational) order. 

In this letter, we consider the role of nematic fluctu- 
ations on the structure of the phase diagram associated 
with the PDW state, and on the properties of its phases. 
Here we focus on the interplay between nematic, PDW 
and the charge 4e SC which we will show also to have 
nematic character. We consider different regimes charac- 



terized by the strength of the coupling between the orien- 
tational (nematic) degrees of freedom of the PDW with 
the underlying lattice, all the way to the decoupled case 
where the system has a continuous rotational invariance. 
In the absence of a coupling to a lattice it is not possi- 
ble to break spontaneously translation invariance in 2D, 
and a 2D continuum smectic is always thermally melted 
by proliferation of (finite energy) dislocations [l6|. Thus, 
the topology of the phase diagram of the PDW state is 
strongly affected by the coupling to the lattice. 

The PDW is an anisotropic quantum liquid crystal 
state that breaks the point group symmetry of the lat- 
tice as well as translation and global gauge invariance. 
We have constructed a phenomenological theory of the 
PDW with a structure reminiscent of the McMillan-de 
Gennes theory of the nematic/smectic transition of clas- 
sical liquid crystals, with several significant differences: 
a) it involves the order parameter of the charge 4e SC and 
its geometric coupling to nematic order, b) the nematic 
order has, in addition to the standard elastic Frank free 
energy [l7[, a coupling that breaks the continuous rota- 
tional invariance down to the (Ising) point group symme- 
try of the lattice. The strength h of the lattice coupling 
defines a temperature T^,, below which, the coupling to 
the lattice breaks the continuous rotational invariance 
down to the point group symmetry (say Z2) of the lat- 
tice, forcing the director to take only one of two perpen- 
dicular orientations. In electronic nematic phases, such 
as those occurring in strongly correlated systems, lattice 
effects are not normally small. Typically nematic order 
occurs on mesoscopic scales, and the effective value of h 
can be smaller than the other couplings. As temperature 
increases, Ising domain walls proliferate and, above Th 
where the domain wall tension vanishes, the Ising order 
disappears and the system becomes invariant under 90° 
rotations. We discuss the behavior and role of the topo- 
logical excitations of this system as h is varied from to 
large values. 

The qualitative structure of the phase diagrams is sum- 
marized in Figsll] [21 and [S] In the regime in which 
Th is large (FigH^), the nematic is the ordered phase 
with higher Tc. In this regime the physics is qualita- 
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FIG. 1. Schematic 2D phase diagrams for high T^: PDW/N- 
SC, Nematic, CDW are KT transitions. The Isotropic- 
Nematic hne is an Ising transition and the Isotropic-CDW 
is first order. 
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FIG. 2. Schematic 2D phase diagram for low Th- Labels are 
the same as in Fig[T] 



tively similar to that of Ref . jl5| . The main effect of ne- 
matic fluctuations is to soften the smectic stiffness (k) 
to Kcff — «/(! + kIQI'^/Zi) where Q is the the smec- 
tic ordering wave vector. As the strength of the cou- 
pling constant h is reduced, the effective CDW stiff- 
ness Keff gets weaker, the PDW portion of the phase 
diagram shrinks and the region of N-SC grows. There 
is also an isotropic/Ising-nematic transition that takes 
place at temperatures higher than the SC transition. 
For large Keff, as the temperature T increases the PDW 
phase may melt (or may not) by a direct PDW to 
isotropic transition or through an intermediate CDW de- 
pending on the sign of the coupling between the order 
parameters in the Landau theory. The PDW to ne- 
matic superconductor (N-SC), CDW and nematic, are 
Kosterlitz-Thouless (KT) transitions, mediated by un- 
binding of double dislocations, vortices or half-vortices- 
single dislocations [ij, [11, [l«] respectively. The direct 
CDW/isotropic transition is likely first order. 

For lower values of Tt we have a softer nematic (Figl2]) . 
Still, Th is an upper bound for the Tc of the PDW since, 
above this temperature, no translational order is possible. 
Similarly, the CDW phase is also absent. 

In the extreme case of a system decoupled from the 
underlying lattice, Th = 0, the otherwise logarithmically 
divergent dislocations have now a finite energy and pro- 
liferate at all temperatures. In this case, the PDW phase 



occurs only at T = as shown in Fig. [3] Provided ps 
remains finite, in this regime the non-SC nematic phase 
also cannot exist since, as we will see below, in the N- 
SC vortices are strongly bound to disclinations, and thus 
proliferate simultaneously. Hence, at Th = the nematic 
phase disappears. In addition, a tetracritical point is 
now possible: nematic and SC order are decoupled and 
the N-SC is a coexistence phase. 
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FIG. 3. Schematic 2D phase diagrams for Th — 0: all lines 
are KT transitions. Labels are the same as in Fig[T] 

a. Order parameters and Landau theory: The PDW 
state has two complex order parameters A±Q{r) ^ 
|A| exp(i^±Q(r)), with Q — Qn and n is the director 
of the nematic order. We now define the SC phase 
field 6 — {0+Q + 9-g)/2 and the smectic phase field 
ip = (e+Q - 9^q)/2 0, [H, [111. In this work, unlike 
Ref. [lal, the director field n = (cos a, sin a) is allowed 
to fluctuate. The coupling of the orientational degree 
of freedom n couples to the lattice through a potential 
V{n) = — (/i/8)cos(4a) with h > 0, that reduces the 
continuous rotational symmetry to a discrete subgroup, 
i.e. Z2 for a square lattice. For h ^ the Ising nematic 
has a critical temperature Th >0, above which the Ising 
degrees of freedom disorder [1^, [20] ■ The structure of the 
phase diagram depends on the values of Th, the super- 
conductor stiffness ps and the CDW stiffness k. 

Aside from standard quadratic and quartic terms in the 
PDW order parameters A±q, the CDW order parame- 
ter pk, and the charge 4e N-SC order parameter A^e, 
the PDW free energy contain trilinear couplings that re- 
late all the order parameters leading, in particular, to 
the relation K = 2Q, and the PDW order induces a 
uniform charge 4e N-SC^^. The phase fluctuations of 
the CDW and the 4eSC are locked to 2(^ and 29 respec- 
tively. In addition, we consider a nematic order parame- 
ter, the symmetric traceless tensor Nij — N('hinj — Sij/2) 
where i,j = 1,2 for {x,y) orientations respectively. For 
T << Th, Nij has essentially two values ±A^ character- 
izing an Ising symmetry. The competition between the 
nematic and the 4e SC is governed by the quartic term 
/3Tr7V^|A4ep. If /3 < 0, the condensation of one phase 
enhances the other. For Th > T^'-^ a multicritical point 
may occur. 

b. Nematic-SC couplings: Local nematic order can 
be regarded as a fluctuating metric gij = dij + XNij de- 



scribing the local anisotropy of the stiffness terms of the 
free energy 
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d = / dx^ v/d^ a'' { (Df^^PK)* (Df^PK) + (1) 



where the covariant derivatives are given by 
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Vi - i2eAi ± iQ 5n., 


(2) 


V. + i2Q 5n, 


(3) 


Vi - lAeAi 


(4) 
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Here, A is the electromagnetic vector potential, and 
(5n = n — no is the director fluctuation about an arbitrary 
direction np. The ± sign in Eq.(l2|) depends on whether 
the covariant derivative acts on A±q respectively. 

The effect of nematicity is, in several aspects, very sim- 
ilar with the effect that a curved surface has on an order 
parameter [21[: due to geometrical frustration, disclina- 
tions can be regarded as representing an intrinsic curva- 
ture of the geometry whose sign is the sign of its topo- 
logical charge of the disclinationJ22j. From Eqs.(IT]), ^ 
and Q, the SC phase field couples to the nematic order 
through the metric, whereas the CDW order parameter 
also couples to the nematic in the standard way through 
the covariant derivative (iTI. [23j . Thus, vortices and discli- 
nations have an attractive ( "gravitational" ) interaction. 

The way that ^n enters in the covariant derivatives 
resembles the McMillan-DeGennes|17l |24| theory for the 
smectic-nematic phase transition in classical liquid crys- 
tals. For small rotations, a « 1, (5n-no = O-t-O(a^) and 
the metric only represents a rigid anisotropy. At this level 
of approximation, smectic and SC degrees of freedom are 
only coupled by topological constraints J15j. This limit 
is suitable to describe smectic-nematic phase transition. 
However, to study the properties of the phase, we need 
to consider the full structure of 5n. Also, note that in 
Eqs. ([3]) and dU, the "gauge fields" couple with the twice 
the value of the smectic and the electric charge. This 
is a direct consequence of global translation and gauge 
invariance that force the CDW wave vector K = 2Q, 
and the superconductor charge to be 4e. The factor 
det(g) = 1 — X^N"^ /A is important near the isotropic- 
nematic phase transition where N fluctuates strongly. 

c. PDW fluctuations: Deep in the PDW phase, we 
take the amplitudes of the order parameters |A±q| = A, 
|A4e| and |pk| = p &s constants. The low energy physics 
is governed by SC {9) and smectic (</?) phase fluctuations. 
Since the temperature T << Th, fio points in the direc- 
tion a = or a = 7r/2. The potential V{hQ + 6n) induces 
a "mass" term ^ /i|(5np. From Eq. ([TJ we obtain. 
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kf + QSni){djip + QSuj) 
+ Ki{V ■ Snf + iCa (V X dnf + h\Sn\^\ 



where ps — A'^ + |A4e[^ is the superfluid stiffness of the 
4e SC, and k = A^ + p^ is the CDW stiffness; Ki and 
K3 are the nematic Frank constants [13|. In what follows 
we will assume that both ps > and k > and that 
they never become small. We can now safely integrate 
over the massive nematic fluctuations. Let no be in the 
X direction. By expanding Sn to leading order in a, we 
obtain in the long wavelength limit 



(6) 



F= I d^x {p, (A, {d^ef + K\dyOf) 



For simplicity, we have set A = 0; Ac and As are (finite) 
anisotropics. The low energy physics is governed by two 
anisotropic XY models where the effective CDW stiffness 
is renormalized to Keff = k/(1 -|- k|QP//i). Hence, at low 
temperatures, the nematic fluctuations soften the CDW 
stiffness and the part of the phase diagram of the 4e N-SC 
is enhanced while that of the PDW shrinks. 

In the limit Th ^ T, Berg et al [I5J derived a phase 
diagram for the thermal melting of the PDW state. In 
this regime, the phase transitions are driven by the KT 
mechanism of unbinding (in this case) double disloca- 
tions, vortices or half vortices bounded to single disloca- 
tions (see also Refs.[l3l [25|). However, for T > Th not 
all of these processes are possible since in this regime the 
free energy of the dislocations of the smectic becomes 
finite and thus always proliferate [ig]. In this case, dislo- 
cations get an energy Ed ~ ln(^/a) where a is the core 
of the defect and ^ is a typical length scale controlled 
by the competition between the Frank constant K3 and 
the CDW stiffness k, ^ ~ \/K3/k. Thus, if Th is low 
enough, all smectic orders (PDW and CDW) are de- 
stroyed and only nematic and the charge 4e N-SC are 
possible. Thus, the PDW state can only exist for T < Th 
and, if Th -> 0, the PDW is suppressed at afl T > 0. 
At T > we have short ranged smectic correlations 
(A+q(x)A;q(0)) - exp(-|x|/0 cos(2Q ■ x), while the 
correlations of the uniform charge 4e N-SC exhibit quasi- 
long-ranged order (A4e(x)A4g(0)) ^ l/|x|''(^), with a 
non-universal temperature-dependent exponent rjiT). 

d. Nematic superconductor: Deep in the 4e N-SC 
phase, above the PDW melting, the free energy becomes 



F= d^x\K\\7a\^ + ps\V0\^- 



^{n-Vef} (7) 



where we have set K^ = Ki = K. The last term in Eq.© 
is due to the coupling between the SC currents and the 
nematic fiuctuations: NijJiJj, where J = ps^O is the SC 
current. This coupling is non-polynomial in the angle a 
of the nematic order, since n = (cos a, sin a). In order 
to minimize the free energy, the current should be lo- 
cally perpendicular to the director n. We find two types 
of topological configurations that minimize Eq.Q: a) 
isolated disclinations, and b) (half) vortices bounded to 
disclinations, in such a way that 3-h = at all points. For 
(5) instance the disclination Ui = Xijr and did = eijXj/r"^ is 



one such configuration. Vortices and discUnations have 
an attractive logarithmic interaction, whose sign is inde- 
pendent of the sign of their topological charges. There- 
fore, once again provided p^ > is that is never small, the 
disordering of N-SC can only be produced in two ways: 
a) by unbinding disclinations, which restores isotropy but 
does not affect the SC (the N-SC/SC transition), or b) by 
the proliferation of (half) vortices tightly bound to discli- 
nations (the N-SC/normal transition). The coupling to 
the lattice anisotropy h changes this scenario. For h > 
the nematic transition becomes Ising like and it is driven 
by the proliferation of domain walls, not by disclinations. 
Hence vortices are no longer bound to disclinations. This 
leads to the phase diagram FigsHEl Notice that a non- 
SC nematic can only occur (in this regime) if ps — > 0. 

e. Quasiparticles of the 4e N-SC: Within a 
Bogoliubov-deGennes approach we expect that a spin- 
singlet bilinear of quasiparticle Fermi fields to couple 
linearly to the pair field A(r). In a conventional SC, 
BCS theory predicts that coupling to lead to a gap in 
the quasiparticle spectrum once the pair field acquires 
an expectation value. However, in the N-SC only the 
charge 4e field has an expectation value. This order 
parameter is equivalent to the condensation of a quartet 
of fermionic quasiparticles. Alternatively, we can view 
the charge 4e N-SC as a state in which a composite 
operator of pair fields has an expectation value, i.e. a 
pairing of pairs. Thus, in the N-SC the pair field remains 
strongly fiuctuating and uncondensed. Since the pair 
field couples directly to the Bogoliubov quasiparticles, 
unlike a conventional (paired) SC, there is no net energy 
gap in the N-SC. Nevertheless, the strongly fiuctuating 
pair field of this phase strong scatters the quasiparticles 
resulting in a reduction of their spectral density, an 



effective "pseudogap" . 

/. External currents: Consider a current J in the 
4e N-SC state. The coupling between nematic and SC 
currents is given by a term of the free energy of the 
form (ps/2) Tr(7VJ), where N is the nematic tensor and 

J'ij — JiJj ~ 2'^^'^y- Thus, a constant current (J), acts 
as an external nematic field that explicitly breaks con- 
tinuous rotational symmetry. This term competes with 
the potential V{n). Integrating over nematic fluctua- 
tions and considering, as before, ha along the x direc- 
tion and expanding Sn{a) to 0{a{x)'^) we obtain the 
eff'ective strength of the symmetry breaking potential, 
/Joff — h + j2{Jy — Jx)- The effect of a current along the 
director is to soften the smectic order. In this way, it 
may be possible to tune the effective CDW stiffness k^s 
and the relative weight of PDW, 4e N-SC and nematic 
phases. 

Here we discussed the phase diagram for PDW/N-SC 
phase transitions as the strength of coupling to the un- 
derlying lattice is varied. So far the only known physical 
system in which a case for PDW order can be made is the 
cuprate La2-a;Ba3;Cu04 near doping 1/8, which would 
put it on the right hand side of FigH] The observation 
of "fluctuating stripe order" and nematic order in under- 
doped cuprates[26] leads us to suspect that a fraction of 
these phase diagrams of FigsfJElmav possibly hide in the 
pseudogap regime, including a form of charge 4e N-SC 
order. These phases may be masked by the effects of 
disorder as "glassy" regimes. 
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